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ABSTRACT 

We present the serendipitous discovery of an L subdwarf, 2MASS J06164006— 6407194, in a search 
of the Two Micron All Sky Survey for T dwarfs. Rs spectrum exhibits features indicative of both a 
cool and metal poor atmosphere including a heavily pressured-broadened K I resonant doublet, Cs I 
and Rb I lines, molecular bands of CaH, TiO, CrH, FeH, and H2O, and enhanced collision induced 
absorption of H2. We assign 2MASS J0616— 6407 a spectral type of sdL5 based on a comparison 
of its red optical spectrum to that of near solar- metallicity L dwarfs. Its high proper motion (// = 
l'.'405±0.008 yr _1 ), large radial velocity (V^ a d = 454±15 km s _1 ), estimated uvw velocities (94, —573, 
125) km s _1 and Galactic orbit with an apogalacticon at ~29 kpc are indicative of membership in the 
outer halo making 2MASS J0616— 6407 the first ultracool member of this population. 
Subject headings: infrared: stars — stars: low-mass, brown dwarfs — subdwarfs — stars: individual 
(2MASS J06164006-6407194) 



1. INTRODUCTION 

Ultracool subdwarfs are metal-poor stars and brown 
dwarfs that have spect ral types later than ~sdM7/esdM7 
dBurgasser et al.ll2005D and T eS < 3000 K (|Leggett et al.l 
2000). Since cool subdwarfs are typically members of the 
thick disk or halo population, they exhibit large space 
motions re lative to the Sun (< V >=— 202 km s _1 ; 
lGizid[j~997l ) and th us are often identified as high proper- 
motion stars (e.g.. iGiclas et aL|[I^7lLll978t lLu vtenlll979L 
iPokornv et al]|2004 ; ILepine k Sharall2005h . Due to their 
intrinsic faintness and paucity in the solar neighbor- 
hood, only ~40 ultraco o l subdwarfs are curren tly known 
([Burgasser et al.1 l2007t ILepine fc Scholzl 120081 ). Such 
small numbers stand in stark contrast to the over one 
thousand near solar-metallicity ultracool dwarfs (spec- 
tral types later than M7 and T cS < 2400 K) that have 
been disc overed in the Two M icron All Sky Survey 
(2MASS; iSkrutskie et all 1200^ . the Sloan Digital Sky 
Survey fSDSS: lYork et al.l l2000f). and the Deep Nea r 
Infrared Southern Sky Survey (jEpchtein et al.l 119971 ). 

Electronic address: mcushing@ifa.hawaii.edu 



many of which po pulate the L and T spectral classes 1 
(|Kirkpatrickll2005D . 

Of the roughly 40 ultracool subdwarfs known, 
only three, 2MASS J0532346+8246465 (hereafter 
2MASS J0532+8246; iBureasser et~all [200l . 2MASS 
J16262034 +3925190 (hereafter 2MAS S J1626+3925; 
Burgasseii[200l . and SDSS 1256-0224 (jSivarani et all 
2009) are classified as L subdwarfs 2 . The spectra of all 



three L subdwarfs exhibit Cs I and Rb I lines, pressure- 
broadened K I lines, CaH, FeH, CrH and H2O bands, 
and enhanced collision induced H 2 absorption (CIA, e.g., 
lLinsky|[l969t iBorvsow et aI1ll997D. all fea t ures indicative 
of very cool (iKirkpatrick et al.l 119991 : iLeggett et al.1 
l2001t iReid et al.ll2001f) and metal-d epleted atmospheres 
(|Besselll Il982t ISaumon et al.l I1994D . In addition, they 
have large proper motions (0'.'6 to 2'.'6 yr _1 ) and radial 

1 A compendi um of known L and T dwarfs is given at 
http: / /DwarfArchives.org 

_f_LSR 1610—0040 was originally classifi ed as an L subdwarf b y 
I Lepine et al.l fl2003T) bu t obs ervations b y Cushing h. Vacca (2006), 
Rciners & Basri (2006), and Dahn et al. (2008b) indicate that it is 
probably a peculiar M dwarf/subdwarf binary. 



2 



Cushing et al. 



"E 

o 

i— 

c: 
o 

"-I — > 

rd 
_c 

~o 

CD 

Q 



+2 



+ 1 







-2 



TABLE 1 
Properties of 2MASS 
J0616-6407 



+2 +1 -1 

Right Ascension (arcmin) 



-2 



Fig. 1. — J-band image of the field centered around 
2MASS J0616-6407 obtained with the CPAPIR on the 1.5-m 
CTIO on 2008-01-17 (UT). The field is 4' on a side and North 
is up and East is to the left. 

velocities (—1 30 to —260 km s" 1 ) c onsistent with halo 
membership. iBurgasser et al.l ()2008f ) recently measured 
the parallax of 2MASS J0532+8246 and confirmed, 
based on its Galactic orbit, that it is indeed a member 
of the (inner) halo. Although no sp ectral classifica- 
tion scheme exists for L subdwarfs, Burgasscr et al.l 
(|2007l l2009t ) have tentatively assigned spectral types 
of sdL3.5, sdL4, and sdL7 to SDSS J1256-0224, 
2MASS J1626+3925, and 2MASS J0532+8246, respec- 
tively, based on a comparison with the red optical 
spectra of the near sola r metallicity L dwarf standards 
(jKirkpatrick et al.lll999h . 

In this paper, we present the discovery of a 
fourth L subdwarf, 2MASS J06164006-6407194 (here- 
after 2MASS J0616-6407) found in a search of the 
2MASS database for T dwarfs. The discovery of 
2MASS J0616— 6407, along with followup near-infrared 
imaging and spectroscopy, and red-optical spectroscopy 
is presented in §2. In §3 we discuss the spectral classifi- 
cation and kinematics of 2MASS J0616-6407. 

2. OBSERVATIONS 

2.1. Source Identification 

2MA SS J0616-6407 was identified as a T dwarf candi- 
date bv lLooper et al.l ((2007) in the course of a search for 
T dwarfs in the 2MASS All Sky Survey Database. It was 
selected based on its very blue near- infrared color (J —K s 
< 0.02) and lack of any optical counterpart within 5" in 
the USNO-A2.0 catalog or the Digital Sky Survey (DSS) 
I and II images. A finder chart for 2MASS J0616-6407 
is given in Figure 1 and pertinent properties are listed in 
Table [TJ Below we describe followup near-infrared imag- 
ing and spectroscopy and red-optical spectroscopy. 

2.2. Near Infrared Spectroscopy 



Parameter 



Value 



Observed 



cv(J2000.0) a 
<5(J2000.0) a 
2MASS J 
2MASS H 
2MASS K a 
Ha (J2000.0) 
Hi (J2000.0) 

Vrad 



6 fe 16 m 40!06 
-64°07'19'/4 
16.403±0.113 mag 
16.275±0.228 mag 
>16.381 b mag 
1405±8 mas yr" 1 
— 51±18 mas yr — 1 
454±15 km s" 1 



Estimated 


d 


57±9 pc 


Vtan 


379±60 km s" 1 


U c 


94±10 km s" 1 


V 


-573±31 km s" 1 


w 


125±53 km s" 1 



a 2MASS coordinates at epoch 
1998.95 (UT). b The read flag is and 
the quality flag is U indicating that 
no flux was detected at the position 
and that the value is an upper 
limit. c We follow the convention that 
the u is positive towards the Galactic 
Center (1=0, 6=0). 

2MASS J0616-6407 was observed with the Ohio State 
InfraRed Imager/Spectrometer (OSIRIS) mounted at the 
4.1 m SOAR Telescope on 2007 March 12 (UT). We used 
the l'.'O wide slit which provides a resolving power of 
B = A/AA ~ 1400 across the 1 to 2.2 /im wavelength 
range in three spectral orders. A series of six 180 sec 
integrations were obtained along the 24"-long slit. The 
AO V star HD 53191 was also observed to correct for tel- 
luric absorption and to flux-calibrate the final spectrum. 
The data were reduced using a modified version (v3.4) 
of th e Spextool data reduction package (Cushing et al. 
2004). Wavelength calibration was accomplished using 
the OH airglow lines in the science frames. The spectra 
from the 3 orders were then corrected for telluric absorp- 
tion using the observations of the AO V st a ndard stars 
and the technique described in IVacca et alj ([2003) . The 
3 spectral orders were then merged into a single spec- 
trum covering the 1.2—2.3 /im wavelength range. The 
final spectrum was then rebinned to one pixel per resolu- 
tion element in order to increase the signal-to-noise ratio 
(S/N) of the spectrum. 

The spectrum of 2MASS J0616-6407 is shown in Fig- 
ure [2] (middle). The S/N of the spectrum is low, rang- 
ing from ^8 in the J and H bands to ~4 in the 
K band, but is nevertheless adequate to confirm that 
2MASS J0616-6407 is not a T dwarf (the J - K s color 
of 2MASS J0616— 6407 corresponds to a spectral type 
of ^T4) since the spectrum lacks the prominent CH4 
band s found in the spectr a of T dwarfs (|Burgasser et al.l 
I2002t iGeballe et aT]|2002D . Also shown in Figure [2] are 
the spectra of the two other L subdwarfs with publ ished 
spectra, 2MASS J053 2+8246 (IBurgasser et al.ll2003f) and 
2MASS J1626+3925 (|Burgasserll2004D . All three spectra 
exhibit blue continua due to enhanced collision induced 
absorption (CIA) of H2 centered at 2.4 fim and broad 
absorption bands of H2O centered at 1.4 and 1.9 /im in- 
dicating that 2MASS J0616-6407 is an L subdwarf. 
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Fig. 2.— Near in fr ared spectra of 2MASS J1626+3925 
(sdL4; IBurgasser etail l2007h. 2MASS J061 6— 6407, and 
2MASS J0532+8246 fsdL7 nBurgasser et ani2003Tl . The spectra 
are normalized to unity at 1.3 /an and offset by constants (dotted 
lines). Prominent atomic and molecular absorption features are 
indicated. The S/N of the 2MASS J0616-6407 spectrum is not 
high enough to detect the 1.25 /im K I doublet. 



2.3. Red Optical Spectroscopy 
2.3.1. Magellan/ 'LDSS-3 

Red optical spectra of 2MASS J0616-6407 were ob- 
tained on 2006 May 10 (UT) with LDSS-3 on Magellan. 
The VPH-red grism (660 lines/mm) with a C/75 wide (4 
pixels) longslit mask was used, with the slit aligned to 
the parallactic angle. This configuration provides 6050- 
10500 A spectra across the entire chip with an average 
resolving power of 7? 1800 and dispersion along the 
chip of -1.2 A/pixel. The OG590 Ion gpass filter was 
used to eliminate second order light shortward of 6000 A. 
Two slow-read exposures of 1800 s each were obtained at 
an airmass of 1.49-1.60. We also observed the G2 V star 
HD 44611 immediately after and at a roughly similar air- 
mass (1.77) for telluric absorption c orrection. The flux 
standard LTT 7987 (a.k.a. GJ 2147; iHamuv et al1fl99l 
was observed on the previous night using an identical slit 
and grism combination. All spectral observations were 
accompanied by HeNeAr arc lamp and flat-field quartz 
lamp exposures for dispersion and pixel response calibra- 
tion. 



LDSS-3 d ata were reduced in the IRAF 3 environment 
(|Todvl[l98l . Raw images were first corrected for am- 
plifier bias voltage, stitched together, and subtracted by 
a median-combined set of slow-read bias frames taken 
during the afternoon. These processed images were then 
divided by a median-combined, bias-subtracted and nor- 
malized set of flat field frames. The LTT 7987 and HD 
44611 spectra were optimally extracted first using the 
APALL task with background subtraction. The spec- 
trum of 2MASS J0616— 6407 was then extracted using 
the G star dispersion trace as a template. Dispersion 
solutions were determined from arc lamp spectra ex- 
tracted using the same dispersion trace; solutions were 
accurate to ~0.08 pixels, or —0.1 A. Flux calibration 
was determined using the tasks STANDARD and SENS- 
FUNC with observations of LTT 7987, adequate over 
the spectral range 6000-10000 A. Corrections to telluric 
2 (6855-6955 A B-band, 7580-7740 A A-band) and 
H 2 (7160-7340 A, 8125-8350 A, 9270-9680 A) absorp- 
tion bands were determined by linearly interpolating over 
these features in the G dwarf spectrum, dividing by the 
uncorrected spectrum, and multiplying the result with 
the spectrum of 2MASS J0616-6407. The two spectra of 
2MASS J0616-6407 were then coadded to improve S/N. 
We measured the radial velocity of 2MASS J0616-6407 
using the Rb I (7948 A), Na I (8183, 8195 A), and Cs I 
(8521 A) lines. The observed line positions determined 
by Gaus sian fits were compared against the rest wave- 
lengths (Ralc henko et al.l 12008, September 191 ) to derive 
a heliocentric radial velocity (after correcting for the 
Earth's orbital velocity) of 480±30 km s -1 , where the 
error in the radial velocities arises primarily from the 
scatter in the line positions. We defer further discussion 
of the spectrum to the next section. 

2.3.2. Gemini/GMOS Spectroscopy 

2MASS J0616-6407 was also o bserved with the 
Gemini Multi-Object Spectrograph (|Hook et al.l l2004j) 
mounted on the Gemini South Telescope on 2007 
September 13 (UT). We used the 0'.'75 wide slit and R400 
grating which yields a spectrum from 5900 to 10100 A 
with a resolving power of R ~ 1300. The OG515 fil- 
ter was used to suppress stray light from shorter wave- 
lengths. A total of five 1800 sec exposures were acquired 
at two different central wavelengths and at two positions 
along the slit. Standard calibration frames including flat- 
field and CuAr lamp images as well as observations of the 
G2 V star HD 60402 for telluric cor rection were also ob- 
tained. The flux standard EG 21 (|Hamuv et all 119941) 
was observed on 2007 September 14 (UT) with the same 
instrument setup. 

The data were first reduced using the IRAF Gemini 
GMOS package. The raw frames were first bias sub- 
tracted (with a bias created using GBIAS) and flat- 
fielded (with a normalized flat field created with GS- 
FLAT). The OH sky lines were then subtracted and 
then spectra were optimally extracted using the GSEX- 
TRACT routine. Wavelength calibration was achieved 
using the CuAr lamp images. The spectra were then 

3 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 3.— Red optical spectra of 2MASS J0616-6407 obtained with LDSS-3/Magellan (black) and GMOS/Gemini (red). The spectra 
are normalized to unity at 8100 A. Prominent atomic and molecular absorption features are indicated. 



flux calibrated using the observations of EG 21 and 
the GSCALIBRATE routine. Telluric correction was 
achieved by linearly interpolating over regions of O2 and 
H2O absorption in the spectrum of HD 60402, divid- 
ing by the raw spectrum of HD 60402, and then mul- 
tiplying the results into the flux-calibrated spectrum of 
2MASS J0616-6407. 

Unfortunately, the resulting spectrum exhibits strong 
noise spikes at the wavelengths of nearly all the OH sky 
lines. Fringing in the detector makes the subtraction 
of the OH sky lines using low order polynomials dif- 
ficult. We therefore re-extracted the raw spectrum of 
2MASS J0616-6407 so that we could first pair subtract 
the images with the 2MASS J0616-6407 at different po- 
sitions along the slit. Although this did not completely 
eliminate the strong residuals, it did improve the spec- 
trum dramatically. Finally we measured a radial velocity 
of +445 ±18 km s _1 using the same technique described 
in the previous section. The radial velocity measure- 
ments from the LDSS-3 and CMOS spectra agree within 
the errors and a weighted average of the two values gives 
454±15 km s" 1 . 

The LDSS-3 and CMOS spectra are shown in Figure 
H The CMOS spectrum has a S/N of 10 to 20 while 
the LDSS-3 has a lower S/N. The spectra exhibit band 
heads of CaH (7035 A), TiO (7053, 8432 A), CrH (8611 
A), FeH (8692, 9896 A), and atomic lines of K I (7665, 
7699 A), Rb I (7800, 7948 A), Na I (8183, 8195 A), and 
Cs I (8521 A). The heavily pressure-broadened K I dou- 
blet, strong CrH and FeH band heads, and Cs I and 
Rb I lines are hallmark spectral features of L dwarfs 
(jKirkpatrick et al.| [l999) while the presence of the CaH 
and TiO band heads are co nsistent with known la te- type 
M and L subdwarf spectra (Burgasscr et al. 2003) ■ Over- 



all the spectra agree well except longward of ^9200 A 
where the LDSS-3 spect rum is systematically hi gher than 
the GMOS spectrum. iBurgasser et al.l (|2007h found a 
similar offset between LDSS-3 and GMOS spectra of ul- 
tracool M subdwarfs and ascribed it to an unknown flux 
calibration error and/or the fact that the 9200 A tel- 
luric H2O band was not corrected in their GMOS data. 
Our GMOS spectrum of 2MASS J0616-6407 has been 
corrected for telluric absorption so we can eliminate this 
as the cause of the mismatch. Since the cause remains 
unknown, we restrict further analysis to A < 9200 A. 

2.4. Near Infrared Imaging 

An image of 2MASS J0616-6407 was taken on 2008 
Jan 17 (UT) using the 1.5m telescope at CTIO with 
the Camera PAnora mique Proche InfraRouge (CPAPIR; 
lArtigau et all [2004) wide field-IR imager in the J band. 
The object was observed close to the meridian at an air- 
mass of 1.2 and with seeing conditions ~1''5. We used 
a 9-point dither pattern with 10" offsets. Each exposure 
was ~60 seconds amounting to a total integration time 
of ^9 minutes. Science frames were first sky-subtracted 
using a sky frame created by median-combining all of the 
science data taken on a given night and then flat-fielded 
using a normalized dome flat. Individual frames were 
shifted and stacked to form the final combined image. 
The reduced science frame was astrometrically calibrated 
using the 2MASS Point Source catalogue. We used a ba- 
sic, six parameter, least-squares, linear transformation 
to obtain an astrometric solution for the CPAPIR im- 
age. The field-of-view of CPAPIR is ^30" so there were 
over 1600 2MASS reference stars to calibrate CPAPIR to 
2MASS. We limited our reference stars to the 450 stars 
with 12 < J < 15 as objects in this intermediate magni- 
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Fig. 4.— Red o pt ical spectra of 2MASS J1626+3925 
(sdL4; IBureasser eTaTI [2007h. 2MASS J 616-6 407, and 
2MASS J0532+8246 (sdL7: IBurgasser et aLl 1200^1 compare d 
to those of 2MASS 1155009 +230706 (L4; IKirkpatrick et al.lll99gt) , 
DENIS-P J1228.2-1 547 (L5: IKirkpatrick eta l. 1999) and DENIS-P 
J0205.4-1159 (L7; Kirkp atrick et al.l I1999Tl respectively. The 
spectra are normalized to unity and offset by constants (dotted 
lines). Prominent atomic and molecular absorption features are 
indicated. 

tude range transformed with the smallest residuals from 
epoch to epoch. The solution reference stars were re- 
quired to transform with total absolute residuals against 
2MASS of < 0.2 pixels (see Faherty et al. 2008 for more 
details). We measure a proper motion of /i Q =+1404 ± 
8 mas yr _1 and l*i a =— 51 ± 18 mas yr _1 over the 9.1 yr 
baseline between the CPAPIR image and the discovery 
2MASS image. The errors were computed by combin- 
ing the residuals of the astrometric solution, which ac- 
counted for 0.16 pixels, with the plate scale of CPAPIR 
(~1" pixel -1 ) under the given baseline. The positional 
uncertainty for 2MASS J0616— 6407 was also calculated 
by comparing the residuals of transforming the X,Y po- 
sition for our target over consecutive dithered images. 
However this uncertainty was negligible compared to the 
contribution from the astrometric solution. 

3. ANALYSIS 

3.1. Spectral Classification 

Determining the spectral type of 2MASS J0616-6407 
is not straightforward because an L subdwarf spectral 
classification scheme has yet to be rigorously defined in 
either the red optical or the near infrared due to the few 
examples known. The depth of the 1.4 H2O band and 
the shape of the H band (see Figure [2]) suggests that 



2MASS J0616— 6407 has a near infrared spectral type 
intermediate between that of 2MASS J1626+3925 and 
2MASS J0532+8246, perhaps sdL4 to sdL6. In order 
to derive a spectr al type based on the re d optical spec- 
trum, we followed IBurgasse r et al.l (|2007l ) and compared 
the GMOS spectrum of 2MASS J0616 -6407 to the spec- 
tra o f field L dwarf spectral standards (|Kirkpatrick et all 
119990 . The best match is the L5 dwarf DENIS-P 
J 1228. 2— 1547 resulting in an optical spectral type of 
sdL5 for 2MASS J0616-6407. Note that a similar com- 
parison in the near infrared is impossible because L 
dwarf spectral standards have yet to be defined at these 
wavelengths (see however Kirkpatrick et al., in prepa- 
ration). Figure Q] shows the red optical spectra of 
both 2MASS J0616-6407 and DENIS -P J1228. 2-1547 
(grey). The strengths of the CrH and FeH band heads, 
width of the K I doublet, and overall 8000-9000 A 
slope of the two spectra agree well. Also shown in 
Figure H are the spectra of 2MASS J0532+8246 and 
2MASS J1626+3925 along with their respective L dwarf 
standards. The strength of the features in the spectrum 
of 2MASS J0616-6407 are intermediate between that of 
2MASS J1626+3925 and 2MASS J0532+8246 in good 
agreement with its derived spectral type of sdL5. 

3.2. Kinematics 

Although the high proper motion and large ra- 
dial velocity of 2MASS J0616-6407 suggest that 
2MASS J0616-6407 is a high velocity star and thus 
possibly a member of the halo population, a dis- 
tance is required in order to compute space veloci- 
ties. We estimated a spectrophotometric distance to 
2MASS J0616-6407 in the following way. First we per- 
formed a weighted linear fit of the 2MASS Mj, M H , and 
Mks values with respect to spectral t ype of the ultracool 
subdwarfs with known parallaxes (IMonet "elaTl fl99l 
Dahn et all2008U ISchilbach et alj|2008t IBurgasser et al l 
20081), LHS 377 (sdM7), LSR 2036+5059 (sdM7.5), 



LSR J1425+7102 (sdM8), SSSPM 1013-1356 (sdM9.5), 
2MASS J1626+3925 (sdL4) and 2MASS J0532+8246 
(sdL7) and find, 



Mj = 8.02+0.313xSpType, Cov = 



M H = 7.77+0.300 xSpType, Cov 



M Ks = 7.44+0.320 xSpType, Cov = 



0.137 
-0.0114 



0.148 
-0.0125 



0.172 
-0.0152 



-0.0114 
0.00106 

(1) 

-0.0125 
0.00120 

(2) 

-0.0152 
0.00151 



(3) 

where SpType=7 for M7, SpType=10 for L0, etc., and 
Cov is the covariance matrix of the fit. For those dwarfs 
with two parallax measurements, we used the weighted 
average of the two values to compute their absolute mag- 
nitudes. Figure [5] shows the absolute magnitudes of 
the subdwar fs along with those of single field M and L 
dwarfs from I Dahn et al.l (|2002f ) and references therein. 
In comparison to the field dwarfs, the L subdwarfs ap- 
pear overluminous in the J band and underluminous in 
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Fig. 5.— The 2MASS Mj, M H , and M Ka values for LHS 
377 (sdM7), LSR 2036+5059 (sdM7.5), LSR J1425+7102 (sdM8), 
SSSPM 1013-1356 (sdM9.5), 2MASS J1626+3925 (sdL4) and 
2MASS J0532+8246 (sdL7) are shown in red. The linear fit to 
these values and the ± 1 a range are shown as solid and dotted 
lines. Shown in black are a sample of field dwarfs from Da hn et al.l 
(2002) and references therein that have not been identified as bi- 
naries in high resolution imaging. All spectral types are based on 
optical spectroscopy. 

the K band. If we assume that the effective temperea- 
ture scal e of the L dwarfs and subdwarfs are similar (see 
howe ver iBurgasser fc Kirkpa trick 2006: Burgasse r et al.l 
we can ascribe this behavior to a change in metal- 
licity. At a fixed T D g and g, L dwarfs become brighter in 
the J band and fainter in the K band with decreasing 
metallicity (|Burrows et al. 2006) due primarily to the in- 
creasing importance of CIA H2 absorption in the K band. 

We estimate a distance of 57±9 pc for 
2MASS J0616-6407 using the absolute J and H 
magnitudes derived using Equations 1 and 2 for a 
spectral type of sdL5 (2MASS J0616-6407 has only 
an upper limit in the K s band). The value is the 
weighted average of the two spectrophotometric dis- 
tance estimates and the error includes the covariance 
terms and an error of ±0.5 a subtype in the spectral 
type. The estimated distance and proper motion of 
2MASS J0616-6407 imply a tangential velocity V tlm 
of 379±60 km s" 1 which is roughly an order of mag- 
nitude greater than the median value o f field L dwarfs 
(jVrba et all l2004t iFahertv et ail I2008D and confirms 



that 2MASS J061 6-6407 is a memb er of the halo 
population. Indeed IFahertv et all (|2008h found no MLT 
dwarfs with Vt an > 170 km s" 1 in their sample of 634 
L and T dwarfs and 456 M dwarfs. We have also com- 
puted (u,v,w) 4 velocities of 2MASS J0616-6407 with 
respect to the Local Standa rd of Rest (LSR) following 
Uohnson fc Soderbloml (|1987l ) updated with the J 2000 
galactic coordinate transformations of iMurravl (|1989h . 
Following the convention that the u is positive towards 
the Galactic Center (/=0, b=0), we find (u,v,w) = (9A, 
—573, 125)±(10, 31, 53) km s _1 after correcting for a 
solar velocity with respect to the LS R of (10.00, 5.25, 
7.17) km s" 1 (jDehnen fc Binnevifl99l . The LSR moves 
at 6lsr=220 km s _1 (0 and y point in the direction o f 
Galactic rotation, 1=90, b=0; iKerr fc Lvnden-Bel]|[l986T) 
in the rest frame of the Galaxy so 2MASS J0616-6407 
has a highly retrograde orbit with 0=— 353 km s -1 . 
The stars in the inner halo of the Galaxy have slightly 
prograde orbits with 0=0 to 50 km s _1 while the stars 
in the outer hal o have retrograde o rbits with 0=— 40 
to -70 km s" 1 (|Carollo et al.l I2007D . The rotation of 
2MASS J0616-6407 around the center of the Galaxy 
clearly indicates that 2MASS J0616— 6407 is a member 
of the outer halo. 

Based on its (u, v, w) velocities we also com- 
put ed the Galactic orbit o f 2MASS J0616-6407 follow- 
ing IBurgasser et al.1 (|2009l ). We examined two simpli- 
fie d versions of the Galactic potential model described 
in iDehnen fc Binnevl (|1998h composed of spherically- 
symmetric halo and bulge mass distributions and 
an axisymmetric, thin exponential disk. The two 
Galactic mass models wer e parameterized according to 
iBinnev fc Tremaind (|2008l Table 2.3), which fit the mea- 
sured rotation curve of the Galaxy but bracket the al- 
lowable ratios of disk to halo mass at the Solar radius. 
The orbit of 2MASS J0616-6407 was integrated using 
a second-order leapfrog method (kick-drift-kick) with a 
constant timestep of 1 kyr over ±1 Gyr centered on the 
present epoch. Energy was conserved to better than 
1 part in 10~ 4 over the full length of the simulation, 
with the error dominated by the resolution of the grid 
on which the disk force and potential were interpolated. 
The Z component of angular momentum was conserved 
to 1 part in 10~ 13 . 

Figure 6 shows the orbit of 2MASS J0616-6407 based 
on the halo-dominant mass model, plotted in the rest 
frame of the Galaxy. The X coordinate is defined to be 
positive towards the Galactic Center to align with our 
defini tion of it, and the Sun is located at (—8-5, 0, 0.027) 
kpc (IKerr fc Lvnden-Bel]lfl986t iChen et al.ll2?)0lh . Pro- 
grade motion in this reference frame is counter-clockwise, 
so the orbit of 2MASS J0616— 6407 is clearly retrograde, 
as well as wide and eccentric, with R extending from 8 
to 36 kpc (e ~ 0.7 from max/min R) over the simulaiton 
period. The orbit of 2MASS J0616-6407 also exhibits a 
fairly broad range of inclinations, with Z spanning ±10 
kpc (max \i\ ~ 15 degrees). For the disk-dominant mass 
model, the orbit is even wider (extending out to 40 kpc) 
but with similar eccentricity and inclination range. The 
breadth of this orbit and its retrograde motion support 

4 We follow the convention that lower case letters refer to ve- 
locities with respect to the LSR while upper case letters r efer to 
velocities with respect to the Sun (Mihalas & Binncv 1981). 
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Fig. 6 . — The Galactic orbit of 2MASS J0616— 6407 over ±1 Gyr, based on the halo-dominant mass model of Binncy & Trcmaine (Model 
I; 2008). The upper panels show the orbit in (X,Y) and (R,Z) coordinates. The X axis is along the line connecting the Galactic Center 
and the Sun while the Z axis denotes height above the Galactic plane. The lower panels show the time evolution of R and Z, with past 
motion denoted by dashed lines and future motion denoted by solid lines. The position of the Sun is denoted with a solid circle. 



the conclusion that 2MASS J0616— 6407 is a member of 
the oute r halo, as this popu lation dominates at R > 14- 
20 kpc (|Carollo et al.ll2007l ). We note that our conclu- 
sions remain unchanged if 2MASS J0616- 6407 is an equal 
magnitude binary. 

4. DISCUSSION 

The Galactic orbit of 2MASS J0616-6407 is dis- 
tinctly different from that of other late-type sub- 
dwarfs (e.g., LSR J0822+1700, 2MASS J0532+8246, 
LSR J1425+7102, SDSS J1256-0224) which all have 
perigalacti cons of <1 kpc a n d apogalacticon s near 
the Sun (fLepine et all l2004t IBurgasser et all 120081 : 
Dahn et alJl2008at IBurgasser et al.ll2009t ). As noted by 



Burgasser et al.l (|2009| ) , these other subdwarfs have kine- 
matic properties which are consistent with membership 
in the inner halo of the Galaxy. 2MASS J0616-6407 is 
therefore the first late- type subdwarf that is a member of 
the outer halo and therefore may represent a new class 
of L subdwarfs. In particular, the inner halo is com- 



prised of stars with a peak metallicity of [ Fe/H] ~ ~l-6 
while the outer halo has a peak near —2.2 (| Car olio et al.l 
120071) . This suggests that 2MASS J0616-6407 may be 
even more metal-poor than the other late- type subdwarfs 
known. However, the techniques needed to derive ac- 
curate metallic ities of late-type subdwarfs are in their 
infancy (e.g., IBurgasser et alj [20071 l200"9l ) and conse- 
quently an accurate measurement of the metallicity of 
2MASS J0616-6407 must await future work. Finally we 
note that the discovery of additional outer halo ultracool 
subdwarfs will be difficult given their intrinsic faintness 
and the large amount of time that they spend away from 
the solar neighborhood. 
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